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In the work the data of the Yakutsk complex EAS array and their comparison with calculation in the case of 
primary nuclei of different chemical elements are presented. The calculation by QGSJET model have been used 
interpreting experimental data. 



1. INTRODUCTION 

It is not possible to measure a mass composi- 
tion of primary cosmic rays (PCR) in the energy 
range 10 15 — 10 20 eV using the direct method. 
Nothing remains, but to resort to the indirect 
methods when for the similar estimates the mea- 
surement of different components of extensive air 
shower (EAS) are used. That can be character- 
istics of longitudinal or lateral development of 
shower in the air. And usually it is connected 
with the analysis of the most sensitive to the com- 
position of shower components which differ from 
each other by the character of formatting and ab- 
sorbing in the atmosphere, for example, of the 
charged particle flux (electron, muon), the flux of 
Cerenkov or ionization radiation. 

2. ENERGY TRANSFERRED TO THE 
ELECTROMAGNETIC EAS COMPO- 
NENT 

Figd presents the EAS Yakutsk array experi- 
mental data and calculations by the model with 
the decelerated and moderate dissipation of the 
energy into the electromagnetic EAS component: 
quasiscaling models (solid line) and QGSJET 
(dashed line) p^. From Fig^Jit is seen both the 
agreement of experimental data and calculations 
by QGSJET model (proton) in the region E > 
3 • 10 18 eV, and disagreement at E < 3 • 10 18 eV. 
The scaling model gives a noticeably greater value 
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Figure 1. A portion of the energy transferred to 
the electromagnetic EAS component by Cerenkov 
light data at the Yakutsk array. 



of E m /Eo in relation to the experimental data 
that is doubtlessly also connected with the break 
of scaling function in the region of ultra-high en- 
ergies. 

The experimental data in Figgis well approx- 
imated by the expression of a form: 
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The relation |Jy is primarily important for the 
comparison of estimations of E obtained at the 
Yakutsk and Fly's Eye arrays 

The calculations in |2] (see Fig^J have been 
carried out by the QGSJET model in the case of 
primary proton and iron nucleus. A good agree- 
ment of our calculations in the case of the pri- 
mary proton is observed. The comparison of ex- 
perimental data with calculations for the proton 
and iron nucleus indicates to the fact that the 
mass composition of particles of cosmic radiation 
in the energy region of 10 17 - 10 18 eV and above 
3 • 10 18 eV must differ. At E > 3 • 10 18 eV the 
mass composition is most likely close to the pro- 
ton one. 

3. A PORTION OF MUONS WITH E th > 
1 GeV 

The showers with Eq > 10 18 eV and zenith an- 
gles 9 < 60° have been chosen. It is required so 
that the shower axis is within the array and not 
less than three muon detectors (one of them is at 
a distance of 1000 m from the shower axis) op- 
erate during the shower. The density of muons 
flux /^(lOOO) is like a median between the den- 
sities adjusted to {R) = 1000 m, according to a 
mean lateral muon distribution function. In this 
case the detectors which have been operated in 
the shower (they have given zero indications), i.e. 
the detector threshold were used in analyze too. 
The muon flux densities at have distances of 300 
and 600 m from the shower axis have been cal- 
culated in analogous way. Results are in Fig[2t 
and 03. Calculations [3] have shown that p^{R) 
weekly depends on a zenith angle in the interval 
of A9 = (0° — 60°). So for the comparison of cal- 
culations and experimental data we take a mean 
angle 9 = 39°. In order to take into account both 
physical fluctuations in measurement of p M (1000) 
and methodical ones (apparatus errors and ac- 
curacy of the locate of the shower axis) we have 
considered the judicial interval in one cr™' 
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. In Fig|2t' andEJs the calculation re- 
sult is shown by dotted line in the case of primary 
proton and by a dashed line for the iron nuclei. 
A solid line shows the upper and lower limits for 
the case if EAS with E > 10 19 eV would gener- 
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Figure 2. a) Portion of the muons with 
Eth > 1 GeV (%). p M (300)/p s (300)(.) and 
p M (600)/p s (600)(o). b) p M (1000) vs. E rela- 
tions for observed events 10 18 — 10 19 eV (square), 
10 19 - 10 20 eV (points) and 10 20 - 10 21 eV (trian- 
gles). Expected ±ler bounds for the distributions 
are indicated for proton, iron and gamma-ray pri- 
mary by different curve as in the legend. 
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ate by a primary 7-quantum. The calculation for 
the primary 7-quantum has been taken from [3]. 
Dots for p M (300) and circles for p p (600) the ex- 
perimental data show in FigEK- I n FigEb the 
showers in the energy range of 10 18 — 10 19 eV are 
shown by squares, the showers with Eq > 10 19 eV 
are shown by dots and the showers of maximum 
energy are shown by triangles. The comparison of 
experimental data presented in Fig. 2a with calcu- 
lations by the QGSJET model carried out for the 
case of the primary proton and iron nuclei con- 
firms the hypothesis on the fact that the consid- 
erable portion of ultimate energy EAS have been 
formed by protons. Their portion decreases below 
the energy 10 18 eV. It is seen from Fig. 2b, the ba- 
sic mass of points (showers with Eq > 10 19 eV) 
falls into the interval for the proton. 23 points 
fall into a zone of superposition of a proton and 
iron nuclei and 19 showers from 116 fall into a 
zone of upper boundary of calculation for the 
primary 7-quantum. It testifies to the fact that 
even taking into account there exists a probabil- 
ity that showers of such energies can be generated 
by neutral particles and, in particular, by a pri- 
mary 7-quantum. Then it is justified to use the 
analysis of directions of arrival of showers with 
Eq > 10 19 eV for the search of the sources of high- 
est energy cosmic rays. In this case the determi- 
nation accuracy of arrival angles of such showers 
must be not worse than (0.5° — 1.5°). 



4. X„ 



FLUCTUATIONS 



A large number of Cerenkov detectors operat- 
ing in the individual EAS events and also the use 
of a new version of the QGSJET model allow to 
obtain quantitative estimations of mass compo- 
sition of PCR. For this aim we have compared 
experimental data (see FigOU and theoretical pre- 
dictions by the QGSJET model for different pri- 
mary nucleus with the use of x 2 criterion. The 
value of \ 2 nas been determined by the equality 
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'expv^max) is the experimental number of 
showers in the AA max interval. N thcoT (X max , A») 
is the analogous number of showers calculated un- 



Figure 3. Depth of maximum distributions for 
difference fixes energy. 
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der the assumption that the mass number of the 3 
nucleus is equal to A i: and P (A{) is the proba- 
bility of the fact that the shower with the energy 4 
E is formed by a primary particle Ai. Then: 

n 

^Vtheor(^max) = ^ P(Ai) ■ -/V t hcor(A max , Ai) (3) 
i=l 

The analysis of form of the experimental dis- 
tribution of A max at optimal value of x 2 with a 
definite portion of probability doesnt contradict 
to the following relationship for five nuclei com- 
ponent: 

E a = 5 • 10 17 eV: p: (39 ± 11)%, a: (31 ± 13)%, 
M: (18 ± 10)%, H: (7 ± 6)%, Fe: (5 ± 4)%; 

Eo = l- 10 18 eV: p: (41 ± 8)%, a: (32 ± 11)%, 
M: (16 ± 9)%, H: (6 ± 4)%, Fe: (5 ± 3)% ; 

E a = 5 • 10 18 eV: p: (60 ± 14)%, a: (21 ± 13)%, 
M: (10 ± 8)%, H: (5 ± 4)%, Fe: (3 ± 3)%. 

Thus in the framework of the QGSJET model 
one can suppose that the mass composition of 
PCR changes transferring from the energy range 
(5 -30) TO 17 eV to energy range (5-30)-10 18 eV. 
At E > 3 • 10 18 eV the primary cosmic radiation 
consists from ^70% protons and helium nuclei, 
a portion of the rest nuclei in the range where 
there is the second irregularity in the energetic 
spectrum type "ankle" doesnt exceed ~ 30%. 
A high content of proton and helium nuclei in 
PCR in the region of formation of "ankle" is most 
likely connected with an appreciable contribution 
into the overall flux of cosmic radiation in the 
Earths vicinity of radiation coming from beyond 
our Galaxy limits. 
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